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Abstract 

We show that the standard model (SM)-like Higgs boson may decay into neutrinos with a sizable 
decay branching ratio in one well-known two Higgs doublet model, so-called neutrinophilic Higgs 
model. This could happen if the mass of the lighter extra neutral Higgs boson is smaller than one 
half of the SM-like Higgs boson mass. The definite prediction of this scenario is that the rate of 
the SM-like Higgs boson decay into diphoton normalized by the SM value is about 0.9. In the case 
that a neutrino is Majorana particle, a displaced vertex of right-handed neutrino decay would be 
additionally observed. This example indicates that a large invisible Higgs boson decay could be 
irrelevant to dark matter. 
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I. INTRODUCTION 


The newly discovered particle at the Large Hadron Collider (LHC) is now identified 
as a Higgs boson [1, 2], Its measured properties such as spin, parity, and couplings are 
consistent with the Higgs boson h in the standard model (SM) of particle physics 3-6] within 
uncertainties, which are not very small yet. Possible deviations from the SM prediction on 
the Higgs boson also have been examined. 

One of those is an invisible Higgs boson decay. Actually, an invisible Higgs boson decay 
occurs even in the SM through an off-shell Z boson Z pair into four neutrinos v, as h —y 
Z*Z* —> 2v2D. Its branching ratio in the SM is of the order of ICC 3 . If once it is found with 
a larger branching ratio than that due to SM processes, this must be a sign of a beyond the 
SM (BSM). Such BSM models include, for instance, a light neutralino in supersymmetric 
models [tJ], a Majoron A. „, ars a. fourth generation neutrino Q and H.ggs porta, 
dark matter jll]. Searches of invisible decays of the Higgs boson h have been carried out 
and to date only an upper bound on the branching ratio of the invisible decay has been 
obtained Q, Hl- 

In this paper, we show that the Higgs boson h would decay into four neutrinos through an 
extra Higgs boson, which can be seen as the invisible decay, in a class of the two Higgs doublet 
model (THDM). The remarkable feature in this scenario is that the invisible final states are 
a SM particle, neutrinos, compared with other BSM models mentioned above where final 
invisible states are new hypothetical particles, such as a supersymmetric particle or dark 


matter. We consider the so-called neutrinophilic THDM |14Ml6], where one Higgs doublet 


provides the mass of the SM fermions, while the other generates neutrino Dirac masses with 
its small vacuum expectation value (VEV). Phenomenology of the charged Higgs boson was 


studied in Refs. 


1?), |l8]. In this paper, we will study a possible phenomenology of neutral 


Higgs bosons in those models. Because of these Yukawa couplings, both extra CP-even and 
extra CP-odd neutral Higgs bosons, H and A, respectively, couple mostly with neutrinos. 
Thus, through interactions between the SM-like Higgs boson h and the extra Higgs bosons 
H(A), as the Z boson makes new decay processes 


h ->• ifWifW oyA^A^ ->• 2u2u 


( 1 ) 


2 




















arise. 1 If the intermediate H or A is off shell, the resultant contribution is comparable to 
the SM contribution by the Z boson and is not so large. However, if either H or A is on 
shell, the resultant invisible decay width is large. 

II. NEUTRINOPHILIC TWO HIGGS DOUBLET MODEL 

The Higgs sector is of the so-called neutrinophilic THDM, where one Higgs doublet $1 
with its VEV V\ generates the mass of the SM fermions, while the other $ 2 generates neutrino 
Dirac masses through its VEV r> 2 •C v\. Such a Yukawa coupling is realized by introducing 
the softly broken Z 2 -parity charge assigned as in Table [I] The Yukawa interaction is given 


TABLE I: The assignment of Z 2 parity and lepton number. 


Fields 

Z 2 parity 

Lepton number 

First Higgs doublet, <Fi 

+ 

0 

Second Higgs doublet, <h 2 

- 

0 

Lepton doublet, L 

+ 

1 

Right-handed neutrino, ur 

- 

1 

Right-handed charged lepton, £r 

+ 

1 

Others 

+ 

0 


by 


Cy — -yi a L a $it Ra — y Ua Q a ®i u R a ~ yd a Q a ^id Ra ~ 2/m A* $ 2 ^ + H.c., (2) 

where $ = icr^*, Q (. L ) is the left-handed 577(2) doublet quark (lepton), and ur, d R , e R , 
and ur are the right-handed (RH) 517(2) singlet fermions, respectively, a denotes flavor 
where we neglect mixing in quarks and i represents the generation index of RH neutrinos. 
If we admit lepton number violation in theory, the lepton number violating Majorana mass 
term 


C M = --u c R .MiU Ri 


( 3 ) 


This possibility was briefly mentioned in Ref. 17] for a very heavy SM-like Higgs boson in a different 


type of neutrinophilic Higgs model. 
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Q. 


also can be introduced jl4J. The scalar potential is given by 






+Ai|$i | 4 + A 2 |d ) 2| 4 + A 3 |$i| 2 |d ) 2| 2 + A 4 |<l ) {d ) 2| 2 + —( < hj < f ) 2 /+ H.c. j , (4) 

where ji \2 is the soft breaking parameter of the Z 2 parity, as introduced above. Conditions 
that the potential^ is bounded from below and a stable vacuum are given by 19 1 


Ai > 0 , A 2 > 0 , 2 a/ A 1 A 2 + A 3 + min[ 0 , A 4 — | A 51 ] > 0 . 

Components in two Higgs doublets, each with a VEV, are parameterized as 


(5) 


$1 = 


S"1 

vi+h\+ia\ 


$2 — 


V2+h2+ia2 


( 6 ) 


s/2 / \ s/2 

Following the concept of neutrinophilic Higgs model, we take Vi ~ v ~ 246 GeV and 
V 2 V\. The smallness of V 2 is due to the small // 2 [14]. We define tan/3 = v\jvi as usual; 
this corresponds to tan (3 1. The states h 1 and I 12 are diagonalized to the mass eigenstates 

(h and H) as 


/ hA 

1 / 

U ) 

1" 


cos a — sm a 
sin a cos a 



(7) 


Because of tan/3 1, <j)f and a\ are mostly eaten by the W and Z bosons, while we can 
identify the physical states as H + ~ A ~ a 2 , H ~ / 12 , and h ~ h\. Then, the mixing 
angle is found to be 


sm a ~ 


V 2 

Vl 


( 8 ) 


Automatically almost, the so-called SM limit sin(/3— a) = 1 is realized. Because of tan /3 1, 

a ~ 0 is realized. The Higgs boson h with the mass of 125 GeV is also composed as h ~ h\. 
From Eq. ([2]), the Yukawa interactions of extra neutral Higgs bosons are written as 


C Y mf Sm °n f H f ~ i— (- cot P)uiAy b Ui - i ^ —cot PfA^f 

' v sm p 11 ' 'i’ 


f=Ui,di/i 


f=di,£i 


— cos a.v a HP R Vi + sin /3u a AP R Vi + H.c. 

a/2 a/2 

We find that H or A decays into mostly neutrinos for 

a/2 mf 


(9) 


Hod 


v tan (3 


( 10 ) 
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Masses of extra Higgs bosons are given by 


m 2 H 

2 . a 3 + a 4 + a 5 2 

— h -2 + 2 v 1 

(11) 

m\ 

2 . A 3 + A 4 — A 5 2 
— ^2 + 2 v ’ 

(12) 

m 2 H ± 

2 . A 3 2 
— ^2 + ~^ V ■ 

(13) 


To be consistent with the electroweak precision test, one neutral Higgs boson mass should 
be close to the charged Higgs boson mass as 

m,H+ — rriA or m H + ~ m h ■ (14) 


Interactions of extra Higgs bosons with h is 


C D — v ( — (A 3 + -^4 — As)A 2 + — (A 3 + A 4 + A5 )H" + h\H + \ 2 


h. 


(15) 


III. EXOTIC SM-LIKE HIGGS BOSON DECAY 


Now we consider a case where either H or A is light enough to be produced on shell by 
the h decay. There are two mass spectra of Higgs bosons that are as consistent with the 
electroweak precision test: 

m H < m h /2 m H + — rriA (16) 


and 


uia < m h /2 -C m H + ~ m H . 


( 17 ) 


From the mass formulas ffTTT) . (IT2]l . and (ll3l) . we find that mass spectra (fTfil) and (fl7jl can be 
realized for 


0 > A 4 ~ A 5 and 0 > A 4 ~ —A 5 , 


(18) 


respectively. 

With couplings (EJ) and (1T5]h if h decays into H or A, which decays into neutrinos, then 
this fraction is measured as its invisible decay. The decay width of h —y HH is given by 


r (h -A HH) = X 2 h 


32nmh 


4 m 2 H 
mi 


1/2 


(19) 
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Br (h -» HH) 



FIG. 1: Contours of the decay branching ratio of h into HH in the A h ~ m-H plane. 


with A H = A 3 + A 4 + A 5 . For the case of h —> AA, we obtain the same result just by replacing 
Aff and m# with A^ = A 3 + A 4 — A 5 and ttia, respectively. The LHC constraints on exotic 
decay modes, h —> HH or h —* AA, indicate that Ah(Aa) < 0(10 2 ) is allowed. We define 


Br(h HH) 


r (h ->■ HH) 
r (h all ) ’ 


( 20 ) 


which is shown in Fig. [T] By combining the constraint on Ah(Aa) and Eq. (TT8|) . we hnd 
that, for the light H (H), 

A 3 ~ —A 4 — (+)A5 (21) 


should be positive and of 0(1), which leads to the deviation in the diphoton decay rate of 
h from the SM value [2(|, as shown in Fig. [2l Here, the brown and magenta shaded regions 
correspond to the region where the extra light Higgs boson is tachyonic and its on shell 
production is kinematically forbidden, respectively. One can hnd that this scenario predicts 
i? 77 ~ 0.9. The signal strength of h —» 77 has been reported as 1.17 ± 0.27 by ATLAS 21] 
and 1.14±g;|| by CMS 22 . 


In the following sections, we discuss more detailed phenomenology which depends on 
neutrino mass nature. 
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FIG. 2: Contours of h — >• 77 rate normalized by the SM value on the A 3 — m^± plane. In the 
magenta shaded region, the lighter extra Higgs boson is too heavy to be produced by the h decay. 
In the brown shaded region, the extra light Higgs boson is tachyonic. Here, the extra Higgs boson 
mass is evaluated with A 4 = — A 5 . 

IV. DIRAC NEUTRINO CASE 


The condition (ITOj) is satisfied by a large tan [3 for the Dirac neutrino case. Thus, the light 
H is, in practice, invisible and we have Br(h —> inv ) = Br(h —> HH) (shown in Fig. [T|). 
The same is true for a light A as well. 

The constraint on the charged Higgs boson, which decays into a lepton and a neutrino, 
is, in fact, stringent. This decay mode is similar to that of a slepton in supersymmetric 
models (l7| and masses of the first and the second generation slepton is constrained as 
m i > 300 GeV by ATLAS 23] or mj > 260 GeV by CMS 24], Although some differences 
due to its decay branching ratio exist |l7j], roughly speaking, there is a similar bound on H + . 
Referring to Fig. [2J we find that a rather large coupling A3 > 3(2) for the ATLAS (CMS) 
bound is required. 
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V. MAJORANA NEUTRINO CASE 


With the presence of the Majorana mass term (j3jl . the neutrino mass matrix is given as 

( 22 ) 


M = 


(loop) y 

ml -J=v 2 


V2 V * 


with the radiative generated mass mu° op ' > 25], 

yakVkpMk 


m 


(loop) _ 


E 


I67r 2 


m 


H 


m H - Ml 


Mb 


In^f 

Ml 


rrr 


m\ - Ml 


In ^ 

Ml 


(23) 


We obtain a light neutrino mass 


Wo? = 


E 


VakVkfP 1/2 


Mb 


(24) 


the mass of a heavier RH-like neutrino — Mb and the left-right mixing angle 9 


26] 


sin 9 ~ 


yv 2 


(25) 


V2M k 

m { l oop) could be indeed dominant-or at least comparable with tree level seesaw mass — 

because of A 5 = ±0(1) from Eqs. (TT8T) and (|2TT) . 

The charged Higgs boson decays into cb or tb, depending on the mass in the neutrinophilic 


Higgs model 


bi- 


The results for t 


le 


section of 1 pb can be found in Ref. 


H± —> tb mode with a normalizing production 


cross 


27]. However, this constraint is not so stringent because 


the actual production cross section is not so large. The LHC data constrains the mass of 
H ± decaying into be between 90 and 150 GeV 28] and H ± decaying into tv 29]. 

We note here one cosmological argument on the Majorana neutrino case. The lepton 
number violation by the Majorana nature of neutrino plays an important role in cosmology. 


Several cosmological discussions on neutrinophilic Higgs model were held in Refs. 


301 


32]. 


One of them is an enhancement A L = 2 washout process by large Yukawa couplings and 
relatively light RH neutrinos in a neutrinophilic Higgs model 30]. Although a discussion of 
baryogenesis is beyond the scope and purpose of this paper, as a necessary condition to have 
nonvanishing baryon asymmetry in our Universe, we roughly evaluate the condition of no 
strong washout of lepton asymmetry, 2 provided a nonvanishing lepton asymmetry has been 


2 This is because lepton asymmetry is a potential source of the baryon asymmetry in our Universe in the 
large class of baryogenesis scenario 
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generated by any means of a higher energy physics process. If this condition were violated, 
it would be difficult to explain nonvanishing baryon asymmetry in onr Universe, because 
any generated lepton asymmetry is washed out. The washout rate is given by Tal =2 - y 4 T. 
The condition Tal =2 < H(T ) at T ~ 100 GeV, with H being the cosmic expansion rate, is 
rewritten as 

y < 10- 4 , ( 26 ) 

which would be regarded as a cosmologically favored region. 3 

Now we discuss the decay of H or A. For m NR < wh/a, an extra neutral Higgs boson 
H(A ) decay produces one light left-handed-like neutrino v and the other heavy RH-like 
neutrino Nr. The amplitude is calculated as 

I M(H/A uN r )\ 2 = 2\y\ 2 (pi ■p 2 ) = \y\ 2 ((pi+P 2 ) 2 -m 2 NR ), (27) 


where p\ and p 2 are outgoing momentum of v and Nr, respectively. In addition, m v is 
neglected and indexes of y are omitted. The decay width is given by 


t(h/a vNr) = \y\ 2 ( m H/A 2 ~ m N R ) 2 - ( 28 ) 

Here, the summation is taken for all kinematically possible modes. An extra Higgs boson 
decays into SM fermions, mostly the bottom quark, through a tiny mixing a. Thus, its 
decay width is strongly suppressed by a large tan f3 as 

< 29) 

Here, we define 


and 


Br(H/A —> inv ) 


T(H/A uN r ) 

T (H/A ->• 66 ) + T(H/A isN r ) ’ 


(30) 


Br(h 2H/2A 2v2N R ) = Br(h HH/AA)Br(H/A -> inv). (31) 


Figure [3] shows the contour plot of the invisible decay branching ratio with thick black lines 
as well as the contour of Eq. (T25|) with thin blue thin lines of H and A, respectively. In 
both cases, the invisible decay branching ratio is large for < 0.1 GeV. The dashed green 


3 One known mechanism of baryogenesis which works without any lepton asymmetry is “baryogenesis via 
neutrino oscillation” 


34, 
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(thick) lines are contours of the typical size of Yukawa coupling y = 10~ 5 (10 -4 ) estimated 
from Eq. (124)1 with the atmospheric neutrino mass. As discussed above, y ~ 1CT 4 would be 
critical when we consider nonvanishing baryon asymmetry in our Universe. In both panels, 
neutrino masses dominantly come from the tree level seesaw at the upper left region, and 
do from m l ° op at the lower right region. In the light A case shown in the right panel of 
Fig. [3j the destructive cancellation between m l ° op and the seesaw term takes place. This 
cancellation makes cuspy curves of the invisible branching ratio and Yukawa couplings in 
contours. For the H decay, = 60 GeV and ttia = 200 GeV are taken. For the A decay, 
m H = 200 GeV and = 60 GeV are taken. 


Br(tf inv) 


Br(A -» inv) 




log(mNR) 


log(mNR) 


FIG. 3: Contours of the invisible decay branching ratio of the lighter extra neutral Higgs boson, 
H or A, for rn v = 0.05 eV on the logj 0 (rriN R /GeV) — log 10 (u 2 /GeV) plane. The thin blue lines are 
contours of sin (9. For the dashed lines, see the text. Here we take mn,^nA = 60 and 200 GeV (left 
panel) and 200 and 60 GeV (right panel). 


Figure 0 is the contour plot of the invisible decay branching ratio of h. Here, = 60 
GeV and rri]y R = 10 GeV are taken. 

The produced RH neutrino Nr decays as Nr —> Z*u, h*v or W*l through a tiny left-right 
mixing of sin# < 0(1O -6 ). Here, a sign of inequality becomes more appropriate as m^ oop ' > 
becomes sizable. For such a left-right mixing of the order of 10~ 6 or less, the displaced 
vertex of Nr decay could be generated, and the decay length of the RH neutrino becomes 
ctn r > 1 cm for rri]y R = 0(10) GeV 36]. For a further lighter rriN R or a much smaller sin#, 
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Br (h -> HH -> 2 v2NR) 



FIG. 4: Contours of the decay branching ratio of the SM-like Higgs boson, Br(h — » 2 H — >• 2n2A%), 
in the A# — U 2 plane. Here, we take mjv fl = 10 GeV. 

Nr would not decay inside the detector. One can see that such a small mixing is realized in 
Fig. 0 

On the other hand, for m/v fi > Wh/a, H or A decays into SM fermions through a tiny 
mixing of a Higgs boson (JHJ) as in the usual type-I THDM. 


VI. SUMMARY 


We have shown that the SM-like Higgs boson could have a sizable invisible decay branch¬ 
ing ratio such as 0(10)% with four neutrinos final states, h —> 2u2u for a Dirac neutrino and 
h —y 2u2 Nr for a Majorana neutrino, in neutrinophilic Higgs doublet models, if one of the 
extra Higgs bosons is light enough to be produced by the SM-like Higgs boson decay. For 
the Majorana neutrino, this becomes a case in the parameter region u 2 1$ 0.1 GeV. Because 


of this mass spectrum of Higgs bosons, the SM normalized decay rate of h -A- 77 is predicted 
to be 0.9. In the Majorana neutrino case, the displaced vertex of a Nr decay also would 
be observed. Although the invisible decay of the Higgs boson was recently widely discussed 


m 


11 ] or applied to dark matter physics 




we emphasis that such a size of invisible 


decay can be realized just within simple THDM without a dark matter candidate. 
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